Joint probability for the Pearcey process 

Mark Adler* Pierre van Moerbeke^ 



Contents 

1 Gaussian Hermitian random matrices coupled in a chain with 
external source 

2 Non-intersecting Brownian motions 



3 The inverse of a tridiagonal matrix and its derivatives [14 



4 Integrable deformations and the Virasoro constraints |24 



5 Integrable deformations and 3-component KP |36 



6 A PDE for the Gaussian matrices coupled in a chain, with 



external potential |37 



*2000 Mathematics Subject Classification. Primary: 60G60, 60G65, 35Q53; secondary: 
60G10, 35Q58. Key words and Phrases: Dyson's Brownian motion, Pearcey process, 
random matrices coupled in a chain, random matrices with external potential, infinite- 
dimensional diffusions. 

Department of Mathematics, Brandeis University, Waltham, Mass 02454, USA. E-mail: 
adler@brandeis.edu. The support of a National Science Foundation grant # DMS-01-00782 
is gratefully acknowledged. 

'''Department of Mathematics, Universite de Louvain, 1348 Louvain-la-Neuve, Belgium 
and Brandeis University, Waltham, Mass 02454, USA. This work was done while PvM was 
a Miller visiting Professor at the University of California, Berkeley, USA E-mail: vanmo- 
erbeke@math.ucl.ac.be and @brandeis.edu. The support of a National Science Foundation 
grant # DMS-04-06287, a European Science Foundation grant (MISGAM), a Marie Curie 
Grant (ENIGMA), FNRS and Francqui Foundation grants is gratefully acknowledged.. 



7 The PDE for the transition probabihty of the Pearcey pro- 



cess 



40 



8 Appendix: evaluation of the integral over the full range |51 



The results in this paper form a step in the direction of understanding 
the behavior of non-intersecting Brownian motions on M (Dyson's Brownian 
motions), when the number of particles tends to oo. Consider n Brownian 
particles leaving from points ai < . . . < and forced to end up at 6i < . . . < 
hq at time t = 1. It is clear that, when n — > oo, the equilibrium measure for 
t ~ has its support on p intervals and for t ~ 1 on g intervals. It is also clear 
that, when t evolves, intervals must merge, must disappear and be created, 
leading to various phase transitions, depending on the respective fraction of 
particles leaving from the points and arriving at the points hj. Therefore 
the region TZ in the space-time strip (x,t) formed by the support (c M) of 
the equilibrium measure as a function of time < t < 1 will typically present 
singularities of different types. 

Near the moments, where a phase transition takes place, one expects to 
find in the limit n/^oo an infinite-dimensional diffusion, a Markov cloud, hav- 
ing some universality properties. Universality here means that the infinite- 
dimensional diffusion is to depend on the type of singularity only. These 
Markov clouds are infinite-dimensional diffusions, which 'in principle' could 
be described by an infinite-dimensional Laplacian with a drift term. We con- 
jecture that each of the Markov clouds obtained in this fashion is related 
to some integrable system, which enables one to derive a non-linear (finite- 
dimensional) PDE, satisfied by the joint probabilities. The purpose of this 
paper is to show, for a simple model leading to a cusp, that the joint proba- 
bilities at different times, do satisfy such a non-linear PDE. The interrelation 
between all such equations and the "initial" and "final" (t ±oo) condi- 
tions are interesting and challenging open problems. Moreover, special cases 
have shown an intimate connection between the integrable system and the 
Riemann-Hilbert problem associated with the singularity. 

The first question is the study of the finite Brownian motion model, 
which, as will be explained in section 2, hinges on the study of the coupled 
Gaussian Hermitian random matrix ensemble 7i„ with external source A, 
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given by coupling terms Ci, . . . , and the diagonal matrix (set n — ki + k2) 

(a \ 



O 



a. 



—a 



O 



V 



—a / 



(0.1) 



The probability of such an ensemble is defined by 

Pn(«, Ci, . . . , Cjn] El, . . ., E^) 

■ ^ J_ 1 ^-^TT(M^+...+Ml^-2ciMlM2-...-2Cm-lMm-lMm-2AMm)^J^^ 

Given a disjoint union of intervals and the associated algebra of differential 
operators 



2r 



i=l 



Ee U[^2?-i,&^?] C M and V,iE,) = ^^ibfr'^, 



l<i<m 



and given the tridiagonal matrix 

/ -1 



J J (ci, . . .,Cm-i) 



i=l 



Cl 



O 



Cl 



-1 



O 



-1 



Cm— 1 



Cm— 1 
-1 / 



(0.2) 
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define the following differential operators: 



A. 



Co 



d 



- T^o{Em) - a- c„_i 



d 



da 



do 



m—l 



-Vo{Ei)+ci^. 

OCi 



Theorem 0.1 The log of the probability P„ (a; Ci, . . . , c^-i; Ei, . . . , E^) sat- 
isfies a fourth- order PDE in a, Ci, . . . , Crn~i and in the endpoints b^^\ feat"* 



of the sets E^, with quartic non-linearity^: 
= (f+CiG- + F-CiGA {F-, F+}c, - ( F+G' 



det 



where 



/ -G+ CiF+ -F+ \ 

G- CiF- -F- 

ClF+ -F+ 

\ CiG- ClF- -F- J 

AfCi log P„ + fc| 1 1 Jim, 

A2C1 ± Jim-^] logP„ T i^/U , 



F-G+)Ci{F-,F+}c-, 
(0.3) 



da 



Ci 



{(C2±2aJi„Ci)^f logP„ , F^} 



At 



with K^i-^ a constant, depending on a, Ci, ki and k2, 



K 



J 



2k^ ijaJmm 



kik2 



a 



Note the robustness of these equations: the equations always have the 
same form (10.31) . regardless of the length of the chain of matrices; only the 
quantities F^, and change, via some minors of the matrix J. 



^in terms of the Wronskians {/, g}x = gXf — fXg. 
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The second question concerns a simple model of non-intersecting Brow- 
nian motions on R and their behavior, when the number of particles tends 
to oo. 

Consider n = 2k non-intersecting Brownian motions on M, all starting at 
the origin, such that the k left paths end up at —a and the k right paths end 
up at +a at time t = 1; see [3 [ISl E]- Inspired by [151 [ISl [7] , the Karlin- 
McGregor formula [T7] enables one to express the transition probability P^" 
in terms of the Gaussian Hermitian random matrices in a chain Pn(a; E) with 
external source, explained above; this will be done in section 2. 

Let now the number n = 2k of particles go to infinity, and let the points 
a and —a go to ±oo. This forces the left k particles to — oo at t = 1 and the 
right k particles to +oo at t = 1. Since the particles all leave from the origin 
at t = 0, it is natural to believe that for small times the equilibrium measure 
(mean density of particles) is supported by one interval, and for times close 
to 1, the equilibrium measure is supported by two intervals. With a precise 
scaling, t = 1/2 is critical in the sense that for t < 1/2, the equilibrium 
measure for the particles is indeed supported by one, and for t > 1/2, by 
two intervals. The heart-shaped region TZ formed by the support of the 
equilibrium measure as a function of time < t < 1 has thus a cusp at 
t = 1/2. The Pearcey process V{s) is now defined as the motion of an infinite 
number of non-intersecting Brownian paths, just around time t = 1/2, with 
a precise scaling; see [71 [HI [TSl E]. The joint probabihty that the Pearcey 
process avoids the windows Ei, . . . , E^ at times Si, . . . ,Sm is defined by 



P 



/p(si)nEi = \ 
\V{sJ n E„ = 0/ 

/ allxj(i±f^) ^ zEi \ 
V alla;,(i±f^) ^^E^ J 



lim P; 

2^0 







where P^° was defined above. The main result of this paper is to show that 
the infinite-dimensional diffusion equation for the Pearcey process can be 
replaced by a finite-dimensional non-linear PDE, which is intimately related 
to the 3-component KP hierarchy and which we now describe. 
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Given Ei := [jZ^[xfti^ xf}] C M, define the space and time gradients 

m 2r^ 



d 



£=1 i=i ^-^i e=i 



_d_ 



space and time Euler operators Xq and Tq and a mixed space-time operator 



m 2r£ 



X. 



d 



e=i i=i 
Theorem 0.2 Then 



dx 



I V{si)nEi = \ 



m 2ri 



d 



[si, ...,s^;Ei,.. .,E^) := logP 



satisfies a 4th order and 3rd degree PDE, which can be written as a single 
Wronskian in the gradient X_x: 



1 



X- 



(0.4) 



+ (A'o + 2To-2) X\(^ - 4(i'_iA'_i - T\)T_iQ \ = 0. 



X- 



In particular Q(s; E) = logP ('P(s) fl = 0) satisfies 

Notice here as well the robustness of the equations. The shape of the equation 
(10.41) is the same, regardless of the number of times one considers. Moreover 
the equations are ^^commutative" : the times and windows can be permuted 
simultaneously. Notice that the term containing X_i is the only one which 
ties up the time Si with the precise set Ei. We expect that this equation can 
be used to derive large-time asymptotics, when t ±oo. Also one expects 
that the PDE's for the sine and Airy processes [1] can be obtained from this 
equation by an appropriate scaling limit. These questions remain challenging 
open problems. 
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1 Gaussian Hermitian random matrices cou- 
pled in a chain with external source 

The present paper studies m Gaussian Hermitian random matrices Mj G Tin, 
coupled in a chain with external source A, given by the diagonal matrix (set 
n = fci + 

(a \ 



A 



O 



—a 



O 



V 



-a / 



lk2 



[1.1) 



and given by the coupling terms Ci, . . . , c^-i; its density is given by 

J_g-iTr(M2+...+M^-2ciMiM2-...-2c^_iM^_iM^-2AM^)^j^^ df^^ (^12) 

For each index 1 < £ < m, consider a disjoint union of intervals :— 
Ui=i[^2i-n ^2?] ^ ^1 ^^'^ define the associated algebra of differential opera- 
tors 



2r 



1=1 

Consider the following probability: 



db 



l<i<m. 



[1.3) 



IP^j (q;, Ci, . . . , Cyti—i ;Ei, . . . , Em) 



m ( Mi-spectrum in 1 ^ 

:= P (Ml, ... , Mm) eflUn: with I ] I 

^=1 M^-spectrum in Em J y 

^ J_ / ^-^Tr{M^+...+M^-2ciMiM2-...-2cm-iMm-iMm-2AMm)^j^^ 

(1.4) 

To be clear, this integral is to be taken over the space of m-uples of Hermitian 
matrices, with Mi-spectrum in . . . , Mm-spectrum in Em , and Z'^ is the 
above integral with all the Ei replaced by M. 



7 



Proposition 1.1 . The following hold^: 

(a; Ci, . . . , Cm-i] El, ... , E^) 
1 



m 

A., n ^^^^ '^'^"^ fj Axtf' 

i=l i=\ 
ki+k2 m 

j=fci+l £=1 



:= iy^-\ . . . , ) and := (^(7],, . . . , y^^^.J. 

The proof of this statement is a standard apphcation of the Harish- 
Chandra-Bessis-Itzykson-Zuber formula, combined with the techniques ex- 
plained in the next section. 



2 Non-intersecting Brownian motions 

Consider n = ki + k2 non-intersecting Brownian motions on M (Dyson's 
Brownian motions), all starting at the origin, such that the k2 left paths end 
up at —a and the ki right paths end up at +a at time t = 1: 



P^" (all Xj{t) e E) 
:= P all Xj{t) G E 



all Xj{0) = 
k2 left paths end up at —a at time t = 1, 
ki right paths end up at +a at time t = 1 



(2.1) 



In the Proposition below we shall be using the Karlin-McGregor formula for 
non- intersecting Brownian motions Xj{t) for < t < 1: 



^Throughout this paper, An{x) = Y[i<i<i<ni^i ~ ^j) ^^'^ Vandermonde 
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/■ 1 

= / ^det(p(t;7i,Xj))i<ij<„det(p(l -t;a;i/,5j/))i<i/j'<„]^o?Xi 

for 

1 (a-i')^ 

p{t,x,y):^^e —. (2.2) 



Consider now the Brownian motions at different times 
= to < ^1 < ^2 < • • • < tm-l 

and set 

Ti = ti+i - ti and — = ^ h - — ^— — , for < j < m. 

Considering m disjoint unions of intervals E£ := IJI=i[^2i-i! ^2?] ^ ^'^'^ 
1 < < m, we show that the two probabihties Pq" and P„, as in (1.4) and 
(2.1), are related by a mere change of variables: 

Proposition 2.1 For = to < < ^2 < ■ ■ • < tm-i < tm < tm+i = 1; 
P^"(a// Xi{ti) e El,..., all Xi{tm) e Em) = Pn (a; Ci, . . . , Cm-l, El,..., Em) 
upon setting 



E,^E,.I , c,^ (2-3) 

[te+i - te)[te - te-i) y {tj+2 - tj){tj+i - tj-i) 



and 



I ^{tm tm-l) tn A\ 



Proof: In the following computation, we shall be using the notation 

^M' _ / M (m)x (m)" _ / (m) (m) n 

X — yX]^ ,...,X^^), X — K-^ki+V ■ ■ ■>■^kl+k2)■ 
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Remembering x, y) is the Brownian transition probability (12.21] . one com- 
putes: 

P±'^(all e ^1, . . . , all Xi{tm) e E^) 

^ L . n 

l<i<n 
1<_7 <m 



^0 ^'^ JEI'x...xi<;^ 



ai,...,afcj^^a 



det , xf^)) det fp(t2 - ti, x['\xf) 

\ ■'J l<i,j<n \ 

... det (p{t^ - tm-i,x^[^'^\ xJ.^M det (p{l - t^, x'[^\ aj 

V ■' / l<i,j<n V 



l<i,j<n V / l<ij<n 

/ TT 7 (?) A / TT 



J "'3-1 b + 



1 /" _ 

ai,...,aftj-^a " J E'^x...xE^ i<i<n j=\ 

1<J <m 

Ifej^+lvi^n— » — a 

det I e '2-«i 1 ... det e det e 

l<ij'<n \ / l<i,j<n \ / l<i,j<n 



n J E^x...y.El^ i<i<n 

1 < J <m 



fcl+fe2 .0)2 ,2.0)J-''+l) 2a.('"' 



(^l!)(^2)!^4' JE'-^x...xE^ 

i=l j=l 

fcl+fc2 _ ^ i2^^^U;^_^a^ rn 

i=fcl+l j=l 
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l<i<n 
l<j <m 



i=l 

ki+k2 



i=fci+i 



by setting, for 1 < i < n, 1 < j < m — 1, 
;i _ yi 



— c,- = -^^ and a = a: 



the change of variables x^^'^ ^— y^^^ induces a change of variables for the 
boundary terms of the integrals 

Ei = EnJ — , 

V (^i 

thus confirming (12. 3p and (12.41) . 

Identity = in the previous set of identities is established by means of the 
following argument, which we explain for the indices j = 1,2: 



2.(2), (3) 

JJ dx[^^ Anix^^^) det ( ) det ( 



l<i<n 
l<j<2 



^2 

l<i,jr'<n \ / l<i,j<n 



n 2.WJ^^ n 2.(2) (3) 



n dx?A^ix^^^^)J2i-irlle^^J2i-irile^ 



l<i<n 
1<3<2 



Upon setting 



1 



•^7r(7r'(j)) -^i ' 
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this expression turns into 



(1) 



Zx . X : 



l<i<n 
l<j<2 



l<i<n 
1<J<2 



Then one uses the symmetry of the integration ranges vis-a-vis i. In general, 
one makes not 2 but m synchronized changes of variables. 
Identity = follows by simultaneously setting 



X 



ki+j 



X 



7r'(j)' 
fci+7r"{i)' 



l<i<ki, l<i<m, 7t' e Sk, 



1 < j < ^2, l<i<m, n" e S, 



k2 1 



subsequently summing over vr' G S^-^, vr" G 8^21 giving rise to A„^(x('")') and 
A„2(x^'"'* ), because of the presence of A„(x^^-'), and then dividing by ki\k2l, 
thus ending the proof of Proposition 2.1. ■ 

When taking the limit for n — > 00, we shall need the following scaling 
(see [7]), assuming k = ki = k2- 

2 1 1 

n = 2k = — , ±a = ±-^7, Xi 1— >• XiZ, tj = -(1 + Sjz'^), for z 0. 

(2.5) 



Proposition 2.2 . Given Ei = IJi=i[^2i-i' ^2?] the following holds: 



F^'^iall Xiih) EEi,..., all Xi{tm) G E„ 



b) = u) 'z 



n = -s 



2 ' 



(2.6) 
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with 



a 



i^j " tj-i){tj+2 - t 



Lj-l 



){tj+2 - tj) 



{1 + SIZ'^){S3-S2) 
(1+S22^){s3-Sl) 



for j 



< for2<j<m-2 



l){Sj+2-Sj + l) 



\/ tT'":J^\~''^M forj = m-l 

\Jtm ^m— 1 2 



6!^) = 6} 



(^) 



2(t, 



2m 

2m 



(1) 



(^) 



■5^ + 1-^*^-1 



/or £ 



2m 



(m) 



7 w T> for 2 < £ < m — 1 

, /or i = m 



Proof: Straightforward from Proposition 2.1, combined with the scahng, 
appearing in (12. 6p . ■ 
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3 The inverse of a tridiagonal matrix and its 
derivatives 

Consider the {k + 1) x {k + 1) tridiagonal matrix, with non-diagonal entries 

Ci, . . . , Cfc! 



/ -1 



J (ci,...,Cfe) 



Cl 



O 



Cl 



Ck 



O 



with 



D{ci, . . . ,Cfc) := < -1 



Ck -1 / 

det J~^(ci, . . . , Ck), for A; > 1 
for /c = 
for k — —1. 



Then one checks, that for 1 < j < m. 



Jlj(ci, . . . , Cm-l) = (-ly ^Ci . . . Cj_i 



D{cj+i, . . . , c. 



m— 1 , 



and 



Define 



Jmjipli • • • 1 Cm— l) ( 1) "'Cj . . . Cm—l 



D{ci, . . .,Cm-l) 
D(ci,...,9_2) 



D{ci, . . .,Cm-l) 



Jlj 



for 1 < j < m 



= for j = 0. 
For later use, we shall need the following identities: 



(3.1) 



(3.2) 



(3.3) 
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Lemma 3.1 

9 JmlJ-mi r) j2 •Jim ^ J n 

ClT^ J — ~^-'ml > C17- — log—— — On , C17- — Jmm — ~^-'ml 

OCi Jii OCi Jii OCi 

d , J ml J mi nil r> T 1 J ml J mi 

Cm-1^ log = 2\-l-2.Jram + 



OCm-1 Jii V '^li 

d , J ml -, 21) (Cj-|-i , . . . , Cm— 2) , c\ Jmm Jii JmlJmi 

log— = -l-^^7 r = ~l~2 

OCm-1 Jii -t^(,Cj+l, • • • , Cm-l) Jii 



Cm—1~^ logt/mj 1 2Ciryi—lJm,m—l ^im 2J.fnm 1 5; 

OCm-1 



im 



Proof: At first notice that 
d_ 
dci 



d 

Ci—D{ci, . . . ,Cj) = -2clD{c3,...,Cj 



2{D{ci,...,Cj)+Dic2,...,Cj)) 



Q 2 

D{Cj, . . . , Cm-l) = ~'^Cm-iD{Cj, . . . , Cm-s) 

OCm-1 

= 2{D{Cj, . . .,Cm-l) + D{Cj, . . .,Cm-2)) 

(3.4) 

and 

i-l 

D{C2, . . . , Cjn-l)D(ci, . . . , Q_2)--D(ci, . . . , Cjn-l)D{c2, . . . , ^-2) = D{ci+i,. . . , C„_i) cl- 



D{ci, . .. , Cm-2)D{Ci+l, ■ ■■ , Cm-l)-D{ci, . . . , Cm-l)D{Ci+i, . .., Cm-2) = D{ci, ■■ Ci-2) . 

(3.5)^ 

Then one checks, by f l3.2l) and fl3.5p . that 



JmlJmi ^i ■ ■ ■ ^m—lJ^(.^lj • • • ; Q 



i-2 



Jii D{Ci+i, . . . , Cm_i)L)(Ci, . . . , Cm-l) 

D{ci, . . . ,Cm-2) -D(Ci+i,... 



D(ci, . . . ,Cm-l) /^(Ci+i, . . . ,Cm_i) 
X _ -D(Q+i, ■ ■ ■ , Cnt-2) 

Jmm j-^r \ • V"'^J 

^[Ci+1, • • • 5 Cm-l) 
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Hence 

. • . , Cm— 2) JmmJli JmlJn 



) Jli 

and by (O, 

Jml ( 1) Cj • • • Cjn— 1 



(3.7) 



Ju D{c 

Moreover, explicit differentiation of (I3.6P and using the identities (I3.5p . (I3.2p . 



one is led to 

d J ml J mi 



Cl 



dci Ju 

-D(Cj+i, . . . , Cm-l)-D^(ci, . . . , 



- -2 V 

and, setting i = m, yields at once the last identity on the first line of the 
statement of Lemma 3.1. Also by fl3.6p and (13.51) 

1 ^ 1 JmlJmi 

-c™_i ^ log 

Djci, ... , Cm-2)D{ci+l, . . . , Cni-l) - D{ci, ■ ■ ■ , Cm-l)D{ci+l, . . . , Cm-2) 
D{ci, . . . , C,n-l)D{ci+i, . . . , Cjn-l) 

107 JmlJmi 

1 ^Jjnm I J 
Jli 

From (13. 7p compute at once 

d , Jim r 

Cl— log-— = dii 

OCi Jli 

and from ([32D, (D and (133|) . 

d , Jim 1 r,-D\'^i+l^ ■ ■ ■ j'^m-l) + D{Ci^i, . . . ,Cm-2) 

Cm-IT. log— = 1-2 — ^ 

OCm-1 Jli J^[Ci+l, . . . , Cm-l) 

,-D(Cj+i, . . . , Cm-2) 



-1-2 
-1-2 



-D(Cj+i, . . . , Cm-l) 
JmmJli JmlJmi 
Jli 
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Also, using 

Z)(ci, . . . , Cm-i) = —D{ci, . . . , Cm-2) — c^_i-D(ci, . . . , Cm-3] 
and the explicit formula (13.21) for J^i, one computes, using (13. 4p . 



d 



Cm— 1 



dCm-1 



log Jmi = Cm-1- log (-1) Q . . . £^-1777 T 

V D{Ci, . . . ,Cm-l) 



Sir. 



dCm-1 

1 2Cm— 1 «-^m,m— 1 ^im 

2 (1 Cjn— 1 t-^mjOi— 1 ) 1 '^im 

/ /^(Ci, . . . , Cm_i) + c^_iD(ci, . . . 
V /^(Ci, . . . ,Cm-l) / 

= 2 -P('^l^ ■ ■ ■ , Cnt-2) _ ]^ _ ^ 

confirming the last formula in the statement of Lemma 3.1. 



Proposition 3.2 For arbitrary z E C, consider the map of Proposition 2.2, 
namely 

((si, . . . , Sm), {u^'\ . . . , u^""^)) ^ ((ci, . . . , Cm-1, a), (&('\ . . . , b^""^)) 
where 



a 



6» 



{l+siz2){s:i-S2) 



i^j-Sj-l)(Sj + 2-Sj + l) 
(Sj + l-Sj_l)(Sj+2-Sj) 



/or J = 1 



for 2 < j < m — 2 



(g^-l-S^-2)(l-S„22) ■ _ ^ _ 2^ 

{Sm-S,n-2){'^-Sm-lz'^) ■' 



[/i(s,2;) 



1 < i < m, 



(3.8) 
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with 



1 {1 + SIZ^){S2 - Si) 

2 V 1 + S2Z^ 



2 V - Si_i 



/or z = 1, 



for 2 < i < m — 1, 



2' 



1 - Sm-lZ^ 



for i = m, 



The inverse map involves the tridiagonal matrix J ^ and can he expressed 
as a fractional linear map in '^"'J^'^'"' for 1 < i < m 



Si{a,c; z) 



1 ^2^4 J^iJ^i _^ 2 

-hi 



f/i(s(a,c; z)]z) 



-az J, 



ml 



-1 



az 



Oi Z JmlJmi ^Jli Jmi \ CX^Z'^ ^\ 2 



A^'ote t/iflt also the entries of J can he expressed as functions of Si: 

Jli 



il-zh;) /(1 + S1Z2)(1 + S2^2) 



^2 - Si 



J m.i 



^z^U,{s, z) 

-{l + z'^Si) I {I - Sm-iz'^){l - SmZ"^) 



^z^Ui{s,z) 
with Ui{s,z) as in 11!^. 8\) . 



(3.9) 



(3.10) 



Proof: Step 1: Upon expanding D(ci, . . . , c^-i) along the j-th row, one 
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checks by and ([S3D, that for 1 < j < m - 1, 

m—l 



n 



_ _ _j_ 

Ci . . . 



D(Ci, . . . , Cj_i)iI'(Cj+i, . . . , Cm-l) — cjD{Cj^2, ■ ■ ■ , Cm-l)D{Ci, . . . , Cj_2) 

D{cj+2, • • • , Cm_i)L)(cj+i, . . . , c„_i) 



m—l 



1 ^(Ci,... 



j D{cj+i, . . . , Cm-i)-D(cj+2, . . . , Cm-i^ 



(3.11) 



n 



1 



and similarly, for 1 < j < m — 2, 



m—l 



rj+2-rj = 



1 

D{ci, Cj)D{Cj+i, Cm-l) - CjCj_^^D{Cj+3, Cm-l)D{ci, 9.2) 

D(Cj+3, . . . , Cm-l)D{Cj+i, . . . , Cm-l) 



m—l 



1 D{ci,... 



-D(Cj+i, . . . , Cm-l)-D(Cj+3, . . . , Cm-l) 



n 



1 

These identities then lead to: 



(rj+i - rj_i)(ri+2 - r^) 

- rj-l)(n+l - n) _ _Jral 

(ri+i - rj_i) J- 



a- , for 1 < i < m — 2 



2 ' 

li 



for 1 < i < m — 1 
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and 

(3.12) 





p2 










1 




J Im 



Step 2: It is easier to show that the inverse map of (13. 9p is given by (13. 8p . So, 
from inverting the fractional hnear map, appearing in (13. 9p . one computes 

and so, inverting this map, one computes, in cascade, 

2 _ 4 1,2 2a^2;^rjJmi 



l-ZSj = — — -, l + zs 



a^zh-iJmi - 2 a^z'^TiJmi - 2 

4a^z'^(rj_i - ri)J^i 

{aHh-iJ^^ - 2)(a22;V,_iJ„i - 2) 

4a2z2(ri_i - ri+i)Jmi 



(a2zVi_iJ„i - 2)(«2^V,+iJ„i + 2)" 
Therefore, using (I3.12p . one checks 

= ^!l^^:^^%^^^ = c?for2<.<m-2 

(1 + siz^){s^ - S2) _ ri(r3 - r2) _ ^2 
(1 + S2^^)(s3 - si) r2(r3-ri) ^ 

(1 '5^^ )('5r7i— 1 ■5^—2) ''"m— 1 ^m—2 ^2 



(1 - Sm-i2;^)(l - Sm^;^) r„-rm_2 "^"^ ^ 

4 Sni— 1 



Z^ (1 - Sm-lZ^)[l - SmZ^) 

and similarly for the expressions Ui in (13. 8p . The signs are all specified from 
y/Si+i - Si and^s^I^T^^. 

Identity (I3.10p is obtained by solving 

Q;2^4 JmlJm» _ 2 
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for iffiiimi.^ by substituting the result in (K9h . i.e., 



Q12: -i^^ 

Ui(s(a, c: z): z) = ——; — 'j^ . 

^ ^ ' ' Q,2 ^4 JrnlJrnA _ 2 

Jli 



and then solving for Finally expressing J^i = leads to the 



second relation (13.101) . The first relation (13.101) is obtained from the ratio 
■Jli = J "^^j^- s-iid using the result previously obtained for J^i at z = 1 

Lemma 3.3 The following identities hold 



a—Si{a, c; z) = -811^ [s, z)— = — ~ s^z 

oa Jml z^ 

d 

a— log Ui{s {a, c; z),z) = -z^Si 
oa 

ci — Si[a, c; z) = 811^ [s, z)J^^ 



d 

ci— log[/j(s(a,c; z),z) = 5^ - 2az^JmiJiiUi{s, z) 

OCi 



^i^-TT^ V(l + + S2Z^){1 - SiZ^) 

2z^{S2-Si) 



Cm— 1 r\ ^i\C>i) C, Z) 

OCm-l 



iUf{s,z)^^^ (l + 2J, 



1-1 —-'mm J 



1 + SiZ 



ml \ ■Jli 
2 



2z (ySffi ■^ni—l) 

((l - SraZ^) (l - Sra^iZ^) (l + SiZ^) - 2 (l - SmSm-lZ'^) (l " SiZ^)) 

d 

Cm-i^ log Ui{si{a, c; z), z) 

OCm-l 

= -1-2— --2a2:^f/i(s, z)J™ l + 2Jmm- 



21 



where 

J mi J ml 



1 -\- ^JiYifi 



Jli 

1 



2z Sm— l)(l ) 

X ((1 - SmZ^){l - Sm-lZ^){l + SiZ^) - 2(1 - SmSm-lZ^){l " SiZ^)) 
-^(Q+I) • • • ) Cm— 2) JmmJli JmlJmi 



(■SjTt Z 1) 

(■^m, "^m— l)(l Z Sj) 
(1 - S„_iz2)(l + s^^2-) 



m— 1 1 



J, 



1 /(I + SiZ^){l + S2;22)(i _ s„_iz2)(i _ s^Z^) 



ml 



JmiJml _ (1 - gm^^) (1 - ^m-l (1 + gj^^) 

Jli 22;2 (s„ - Sra-l) (1 " gi^;^) 

-{^ + s,z^) 



J mi Ui 01 
JliUi 



2^3 

l-SiZ^ \{l + SiZ^){l + S2Z^) 



Az^ V S2 — si 



Proof: Differentiating the first identity (13.91) with regard to a and using the 
second expression fl3.9p yield 

f \ Jli QTT2f \-JmiJli 

a—Si{a, c; z) = — = -SU^ (s, z)— 



Differentiating f/j as in ( 13.90 . and using the expression ( 13. 9p for Sj, yield: 

Ot— log f/i(a, c: z) = = — z Sj, 
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while f l3.9p and Lemma 3.1 yield 



d . , 4o Z d JmlJ-mi 

—Si{a,c;z) = T^ci- 



2 



8Ui{s,zfj'^ 



li 



Cm— 1 o Si [CX^ C, Zj 2 Cm— 1 o j- 

CCm-l ( ^2 JmiJmi ^4 2 ' t^Cm-l -/li 

,7 7-/ \2'^miJli I I o T J ml J mi 

iUi{s, Z) — I 1 + 2Jram J 

•-'ml \ '-'li 



d d Jml a^Z^ d JmlJmi 

c,g^logC7.(.;.) = c,-log— --53_-^c.^-^ 

= 6ii - 2az^JmiJiiUi{s,z) 



d 

Cm-i- logf/i(s; z) 

CCm-1 

^ ^ lo '^'"^ Oi Z ^ (9 J ml J mi 

"""^SCm-l -/li Q,2 Jml J,n. ^4 _ 2 "~^gCm-l ^li 

_ ^ ^ -P(Ct+l, • • • , Cm-2) 
D(Ci+i, . . . ,Cm-l) 



"2Q!Z Uii^Sj Z^Jmi ( '^Jmm 



J 



li 



yielding most of the differential identities of Lemma 3.3. The remaining 
relations are a consequence of (13.81) . (13.91) . (I3.10p and (I3.13p . ■ 
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4 Integrable deformations and the Virasoro 
constraints 

In order to compute the differential equation for 



P„ (a; Ci, . . . , c^— 1 ;EiX ...X E^) 

^ P ki+k2 m 



ki 



1=1 
ki+k2 



i=fei+l 

we need to add to tlie numerator of P„ many auxiliary variables 
t := (ii, t2,...),s := (si, S2, ■ ■ ■)> ^ — ^2, ■ ■ ■) and 
:= (^f ), . . .) for 2 < £ < m - 1 
. := (cS)p,g>i for 1 < £ < m - 1 



(4.1) 
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Note that the time variables t, s, u are totally different from the t-variables 
appearing in the Brownian motion, yielding the following integral^ 

n,k, (i, s, u- P, 7(2), . . . , c(i), . . . , c(™-i), a,E,x...xErn) 

„ ki+k2 / m 



fci 



j=l 



(m) (m) 



A:i+A:2 

^ ' '''tj)l<i<ki, l<i<fci+fc2 

det I I (4.2) 



if^ij)l<i<k2, l<j<ki+k2 



where 



m 



_ 1. 

e 2 



i Er=i x(*)2±ax(™)±/3x('»)2+Ep,,>i E.T/ 4?(xW)''(x(^+i))«+Er="2' E.°ii 4'\--^^'^r 



(4.3) 

This is to say, the integral (14. 2p . along the locus 



C = I *i = 0' *i = 0' = 0' = 0' T^^^ = 0' ^ (4 4) 

^ cfl = Q and cjf = for i,j>l with (ij) ^ (1, 1) ' ^ ' ' 



■^If m = 1 or 2, the formulae below must be reinterpreted; e.g., for m — 1, the c^'s are 
absent and for m = 2, 7 is not present. 
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yields the integral (14. ip : also for the sake of brevity, set 7^ := '-)f\ The 
following locus 



ti = 0, Si = 0, Ui = 0, 7, 



0, 



Cn = Q and c^-^'' = for i,j>l with (z, j) 7^ (1, 1) 



(4.5) 



will also be used. Then the following statement holds: 



Proposition 4.1 Given a disjoint union of intervals and the associated al- 
gebra of differential operators 



2r 



i=l 1=1 

the integrals satisfy, besides the trivial relations 



-1 d 



db 



{£)■ 



d ^ d d 
dsi dui (9a ' 



d d d 



and, upon setting 



7i := ti, and 



d 



ds2 du2 d(3 



d d 



djm ' dsi dui 
the following Virasoro relations (2 < i < m — 1): 



f 9 

d~ii 



k2 



oti d72 



d 



ic. 



(1) 



d 



i>2 * ^ i>2 

J>1 



8r^ ' 



E 

i>2 



1< i<m 



(4.6) 



.(1) 9 

'^^i a (2) 
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d d ^ /, , N 
Ci-i— — + V [h + hnt 



e-1) 



d 



d 



t>l L/Lj j_l i>2 OfC-j-l 

J>2 J>1 



+ Er>2 n-^ 



w d 



lfe2 




—ki{si — a) — k2{ui + a) 



i>2 



d 



d 



Ui 



Jm-l) ^ 



Q- (m-1) 



. (m-l) 

o^(m-l) 



i>l 
J>2 



9c, 



(4.7) 



an(i ( only needed for i — l,m) 



^ (/a + A-2)(A-i + A-2 + l) \ 
'5ci 2 



V 



(<.3)7^(l.l) 
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da. 

, d d \ d , d d . 



dcm-i ds2 du2' 



90 



^ ds2 0U2 2 2 



d 



d 



9 



(m-l) 



(4.8) 

Before giving the proof of Proposition 4.1, we need the following lemma, 
concerning the expressions: 



din := A„(x) JJ^ dxkC^i 

k=l 

n 



Lemma 4.2 



fe=i 



0, k = -l 



5 



e=0 



d I n(n+l) 



dln{x), k — 
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J>2 



iti 



e=0 

d 

dii-i 



d 



+ nti 



i>2,i>l 

d 



ds 



k = 



,i>2 



e=0 

d 

dsi-i 



d 



+ nsi+ -^'^^i^^^^ + E^iiJ" ) ^(^'^)' ^ = -1 



«>lj>2 
5 



i>l 



A; = 



Proof: This is obtained by setting i— >• Xj + ex^'*'^ and then yi ^ yi + eyf^^', 
then take the derivative with regard to e and set e = 0. 

Proof of Proposition 4-1- 

Case 1 : Performing the infinitesimal change of variables, for all 1 < i < ra, 



(1) (1) , (1) 

X ■ I — ^ XI + exi 



fc+i 



in the integral fl4.2p involves the following integral only 



fcl+fc2 



i=l 



1=1 



1=1 



and leads to the first Virasoro constraints in (14. 7p . upon differentiation with 
regard to e, setting e = 0, applying Lemma 4.2 and taking into account the 
variation of the boundary term in the integral, ll^p ^ b^p —e{bj^^)^^'^-\-0{e'^). 
The shifts in the time parameter produces the terms in the first line of the 
Virasoro constraint. 
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Case 2 : Performing the infinitesimal change of variables, for all 1 < i < n, 

xl I — > xl + exl 
in the integral (14 ■2p involves the following integral only 

i=l 1=1 

and leads to the Virasoro constraints for 2<£<m — lin ( 14. 7p . upon 
differentiation with regard to e, setting e = 0, applying Lemma 4.2 and 
taking into account the variation of the boundary term in the integral, b^j'^ i— >■ 

b^^^ — e{b^pY^^ + 0{e'^). Also here one must take into account the shifts. 
Case 3 : Performing the infinitesimal change of variables, for all 1 < i < n, 

(m) (m) , (m)^+^ 

XI I — > x\ + exl 
in the integral (14. 2 p involves the following integral only 

fcl+fc2 

JJ^ dx^ Ak^{Xi , . . . , ) Ak2{^ ki+l-: ■ ■ ■ ^ ki+k2) 

i=l 



i=l 
ki+k2 



i=ki 



and leads to the last Virasoro constraint for i = m with a similar argument; 
this ends the proof of Proposition 4.1. The relations (14. 6 p follow at once by 
inspection of the integral (14.20 . ■ 
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Setting 



/ -1 



J ^ J ^(ci, . . . , Cm-l] 



Cl 



o 



Cl 



O 



V 



Cm— 1 



Cm— 1 
-1 / 



define the differential operators Af,Ci and .4i,Cj: 

.7 = 1 ^ ^ 



(4.9) 



:^^ + /3J, ^ 



9q; 



- J]JiJ^^-i(&^''^)-25,m/9^J =:Ci-2/?Ji^ 

.7=1 ^ ^ 



da 



9 



Cm— 1 



5Cm-l 



C2 := 
We now state: 



d 



OCi 



Proposition 4.3 Along the locus C, the partials and second partials of f : = 
logTfc^fcj with regard to ii,Si,Ui can be expressed in terms of the operators 
Ai,Ci and 8/8(3: 
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Atf 
Cif 



A2 + 



ld_\ 



2d (3 



f 



AtCJ 



df a (^2 - h) J 

OSi 2 

df a {k2 - ki) 

~r 'J mm 

OUi 2 

- a(fc2 - /Cl) Jml 

OTi 

|^ + 7(^:2' + fc2 + A;i2 + fci) 
<9s2 4 

^ + ^(fc2^ + fc2 + fcl^ + fcl) 



C2/ = ^-^(A;2 + fci)(A;2 + A;i + l) 



dtidsi 
dtids2 



k-lJlm 
- a{k2 - ki)JimJn 



ar a- ~ "(^2 - ki)j'^ 
dt2dsi 



Im 



(4.10) 



Proof: The Virasoro relations (14.71) can be written as (remembering ti = 71) 



P_i(6«) 



V V_,ib(-^))-2p£-ik^-k2)a J 



f = J-' 



/ ^ \ 



f + 0{Cp), 



remembering the locus as in (14. 5p . and the tridiagonal matrix J ^ de- 
fined in (14. 9p . The symbol 0{Cp) means a term which vanishes along 
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Therefore, using the trivial relations (14. 6p . one finds 



J 



/ 



which combined with 

d d 

+ 



a 

d 
972 



a 



a-jm-i 
a a 



\ 1 / 

\ (9si du\ / 



d 



dsi dui da ' 
leads to the following equations for / = logrfc^^j along the locus Cp: 

Of a 

•^l f = TTl 1- 77 (^2 — ki) Jmm 

OUi 2 

df 

Clf = a(^2 - ki)Jml 

Oti 

The second set of Virasoro relations, combined with 

d d d 
ds2 du2 d/3 ' 



leads to the following equations, also along the locus C 



/3) 



df 1 

Atf = ^ + -{k2^ + k2 + ki^ + k,) 

OS2 4 



C2f 



df 



+ -{k2^ + k2 + ki^ + ki 



du2 4 
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Moreover 

= „- „_ kiJim — PJmmJml{k2 — ki) 

OtiOSi 



Indeed, one first needs to check 

1 



Ai{a) = JmmP 2' 



A2{aJmi) = -- a— + Cm_i- ] aJmi 

2 \ da dc^^ij 

J ml J mm 



C'2^Jmrn) '-'1 o J mm ^'-^ml' 



and 

d_ 

dci 

Then, using the expressions for the Aff and Cif, one computes along the 

locuci jCf3, 



AfCJ = At{^^-a{k2-k{)J^^ 



^ytf)-{k2-k,)Jmi{Ata) 

7^ TZ- + 77(^2 — kijJmm — —Jml{kl + ^^2) — (^2 " ki)Jml{JmmP 
Oti \OSi 2 2 

Jmlkl — PJmmJml{k2 — A^l) 



dtids 



^Since one takes the derivative with regard to ti in the expression below, one must 
keep track of the term containing ii in the identity for f before setting ti — 0. 
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and 



f df 

A2Cif ^ ( ^ - «(^2 - kl)Jml 

d 

= g^A2f - A2{a{k2- ki)Jmi) 

— q:^ q- f — (^2 — ki)aJmlJmm 



C2Atf = ^C2f + C2{-{k2-h)J^^) 

osi 2 



dsidt^ 



a{k2 - J^i- 



Finally, set = in each of these expressions. That means that we can 
replace all Ai, Ci by Ai^ Q, except for A2 and B2 , which gives an extra 
d/dp and the composition 



/3=0 



This estabhshes Proposition 4.3. ■ 

Remark: The five operators Af, A2, Ci, C2 form a Lie algebra (upon using 
Lemma 3.1): 

[A^A] = 0, [At, Al] = 0, [A, C2] = 
= -Jim{At + A-^), [C2A] = -(l + 2Jii)(Ji 
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5 Integrable deformations and 3-component 
KP 



The following integral 

n,k,{i, s,u-a, P, 7(2), . . . , 1), c(i), . . . , c^^-'\ E,x...xE^) 

({t^ij)l<i<ki, l<j<ki+k2 \ 
(5.1) 
{fJ'ij)l<i<k2, l<j<ki+k2 ) 

where 



^3-i^j:Th^A (5.2) 



with regard to the inner-product (m > 2) 

Jute. 



for 



1 

From (|1]), it follows that the function Tk^k2 above, which is expressed as 
the determinant of a moment matrix with regard to two different weights, 
satisfies the 3-component KP and thus it satisfies in particular, the following 
PDE's: 

, Tkx+lM'^k^-lM 
or ^Og^fcifc2 — 2 

, Tfci,fc2+l'''fel,fc2-l 

a- ^Qg^fcifc2 = 2 



, _ ''"fci+l,A;2-lTA;i-l,fc2+l 

'k^M 



, • KH-i,K2 — i ' Kl-i,«.2i-i /r- o\ 

A- Q— l0g^fclfc2 = 2 

OSidUi Tf 
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and 



'kik2 



d 



'O^^ki+lM j Tki-l,k2 ~ Tki+IM \ -Q^'Tki-l,k2 



d 



92 



ds2dti 
92 



kik2 - 



d 



-Q — Tk^-IM j ^fci+l,fc2 -'^ki-lM \-g—Tki+l,k2 



d 



ds2dui 



— logTk,k2 



'fclfc2 



d 



^^Tfci-l,fc2+l I Tki+l,k2-l " ''"fci-l,fc2fl I ^^''"fcl+l,fe2-l 



d 



(5.4) 



Thus, upon taking the ratio of the first, second and third equations of (15.41) 
and (15. 3p . one finds 



a 1 



Alog^^iii^ 
<9si Tfc^-i^fca 



ail 052 ^°S'^fci:'=2 
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dtids 



-\0gTk,,k2 



d 1 

^ ''"A:i+l,fe2-l _ duids2 ^^S^k-tM 



dSi Tk^^ 



■l,fc2 + l 



92 

(9ni9si 



logTfc, 



(5.5) 



fc2 



6 A PDE for the Gaussian matrices coupled 
in a chain, with external potential 

The purpose of this section is to prove Theorem 0.2. Notice the probabihty 



P„ (a; ci, . . . , Cm-i] Eix ...X E„ 
1 



l+fe2 



ki 



1=1 



i=l 
ki+k2 

A.,(x('")") n e 



(6.1) 
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is invariant under the involution 

L : a < > —a, (3 < > — /3, Ui < > Sj, ki < >• 

Proof of Theorem 0.1: The first equahty in each of the expressions (16.21) 
and (16.31) below follows from the expressions for Ai f and Cif in (I4.10p . 
whereas the second equalities follow from (15. 5p and the third equality from 
the expressions for A^Cif , ^^Ci/, C2A1 f in (14.101) : 

Ci log = log h 2aJirr, 

— g2 — : h ^aJir 



a^logrfc,,fc2 



2 



92 1 



(6.2) 



log = — log aJ^rn 

Tfci-l,fc2 Tfc^-l,fc2 



92 1 



H+ -X 



F+ 



- aJmm (6.3) 



where F"*", X, if^^ can are functions of logr^j^j, which can also be expressed 
in terms of the actual probability P„, taking into account Lemma 8.1, 

^+ - logrfcj,fc2 = A^Ci logrfcj,fc2 + kiJim = A'^Ci logP„ + fci Ji„ 



dtidsi 
1 d 1 d 

Ht = ^l°g-'=.-'=2 + ^^C,logr.,,,, 

d\ 

- d \ kik2 

A2C1 + JlmT;~ logPn — 2akiJimJmm H Jim 

oa I a 
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= - (C2 + 2aJimCi) At log Tk^^k2 - oi{ki + k2) {Jlmf 
= -(C2 + 2«Ji^Ci)^+logP„ 

Thus, Ai of the right hand side of (16.21) must equal Ci of the right hand side 
of fl6.3p and furthermore noticing that CiaJmm = 0, one finds: 



or in Wronskian notation 

[X.F^]^^ = {Ht,F^}c, - {Ht,F+U = (6.4) 

and 

- {^'^"}c-. = {Hi^F-]c. - {H2.F-]^^ = G- (6.5) 
upon applying the involution l to the first equation. Equations (16.41) and 
(16.51) yield a linear system of equations in X and CiX, leading to 

G-F+ + G+F- 



X 
CiX 



-F~{CiF+) + F+{CiF- 
G-(CiF+) + G+(CiF' 
-F-{CiF+) + F+{CiF~ 



Subtracting the second equation from Ci of the first equation yields the fol- 
lowing: 

[f+CiG- + [f+CiF- - F-CiF+^ 

-(^F+G- + (f+CIF- - F-ClF+^ = 0. 

The second way of expressing these equations is to write a system of 4 
equations, consisting of the system (16.41) and (16. 5p and that same system 
acted upon by Ci. 

= + 
= -{X,F-}^^ - G- 
= -Ci{X,F+}^^+CiG+ 
= -Ci{X,F-}^^-CiG-, 
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or, in matrix notation. 



/ 


G+ 


CiF+ 


-F+ 







1 1 \ 




-G- 


C^F- 


-F- 







X 




CiG+ 


CfF+ 





-F+ 




CiX 


\ 


-CiQ- 


ClF- 





-F- 


) 





= 



and thus the matrix must be singular, leading to the second formulation and 
ending the proof of Theorem 0.1. ■ 



7 The PDE for the transition probabihty of 
the Pearcey process 

Given Ee := ULil^S-i' 4?] ^ M, define the " space'^ and ''time'' operators 
Afc and related to Brownian motion, together with a mixed "space-timd' 
operator ALi, 

e i=i ox^ £ * i i=i ox^ 

We shall also need the intermediate operators 



i=i dxf 

Define a new function 



logP„(a,ci,...,c„_i;6«,...,6("^)) 



(1) ^M^ 



2 



— ■ ■ ■ ) , . . . , X 



(7.1) 
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by means of the change of variables, defined earher 



a — - 



for j 



(l+^l2-^)(g3-^2) 
(l+S2Z^)(s3-Sl) 



(^i-^j-i)(^j+2-^j+i) 

(Sj+l-Sj_l)(s3+2-Sj) 



for 2 < J < m — 2 



(g,„-i-s^_2)(i-^mz^) for 7 = m — 1 



' 2x(^) 



2x, 



(^) 



2x) 



(m) 



(1 + ,S1Z2)(,S2-Sl)' 

''■'i + 1 ~ -'' < - 1 
(Sf-S£_i)(s£ + i-S^) • 



(Sm Sm— l)(l Sm^*^) ' 



for £ = 1 

for 2 < £ < m - 1 
for £ = m 



with inverse map 

1 a^/isii-iiii + 2 



with 



Then 



^^k = bk'Ui{s{a,c-z)-z), 



Ui{s{a,c]z);z) = — 



1 / az ^^hll^ 



mi \ Ct: ^ — 



^ _ ^(^fc^^») _ _ Ai)rj d\ogUi _ (i) dlogUi 

da da ' da ^ ' da ' da 

dxf ^ djb'i^Ui) ^ ^{{)dUi ^ ^(^)^^ d log Uj ^ ^jj^ dlogUj 
dcj dcj ^ dcj ^ * dcj ^ dcj 
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Thus, setting fl7.3p in the right hand side of fl7.ll) and using the chain rule, 
one computes 



(91ogP^ 



da 



da dsi ^ ^ da f)^^^ 

* i=l fc = l LZ-^A; 



i=l 

X 



III' c~\ r^ III' r^ 1 7-7- ^' i r^ 

da ds, ^ da ^ 

j=i * 1=1 k=i ^-^fc 

^ fdsi d dlogUj \ 

1=1 ^ ' 



dlogP„ 

dCj n= 



2 



dcj ds 

1=1 



i=l fc= 



^ dco 

1 ^ '-^■^k 



(i) 



E™" dsj d ^ d log t/, 

j=i 



2r 



dc,- ds 

, j=i 



(0 



d 



J fc=l LZ-^A 



^ dc,- ds 
1=1 ^ 



P_i(E,)logP, 



dlogP„ 



2 

= 7T 



d6' 



2ri 

E 

fc=i 

k=i '-'"k '^•^k 

2r,. g 



2 

= 7T 



fc=i f^a;[.^ 
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2r,: 



Do (^i) log P. 



k=l 

= 



96; 

Si) 



(i)91ogP„ 
k=i ^"k ^-^k 



2 



2ri 



= E 



X 



d 



k.. ' 



= A'o(x«)Q,. 
Using the information above, one computes {e — ±1) 

^flogP„(a,Ci,...,c„_i;6«,...,6(™)) 



n=-% 



1 ( ^ d\ 
1 / 



_ 2 

"77 



9s,- (9 (9 log U4 



da ds-i 



da 



Aid 



(7.4) 



CilogP„ 



^JyP_i(Ej)logP„ 



(7.5) 
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^logP. 



n=4 



da J dsj 



I 



V 



« — 1- OL- 

dlogU. 



CslogP, 



aiogc/. 



E/ 5s, d 



9ci 9s,- 



dci 



Lemma 7.1 The following expansions hold near 2; ~ 0; 



a 



dsj 

da 

dsj 

da 



^ 22 

^2 



d 

a—\ogUi{s{a,c;z),z) = -z^Sj 
— \ogUi{s{a,c]z),z) = 



Cm— 1 



(9 

Cl-R — Si 
OCi 

d 

^ 

OCm-l 



2(S2 - Si) \z 



Sjji—i \Z 
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1 + ^1 + ^2^-2^^ +0(1) 
z 

1 ~l~ ^m— 1 2Sj ^ 0(1) 



2;^ 



d 1 



J-miUi — i ^ + —{Sm + Sm-l — 2Si) + 0(2:^^ 

JuUi = . , ( \ - + S2 - 2Si) 

o 

2 ^2 

a = -y/Sm - Sm~l (1 + — (Sm + Sm-l) + 0{z^)). 

z 2 

Proof: These series follow from expanding the expressions given in Lemma 
3.3 near z ~ 0. ■ 



Lemma 7.2 The operators A\, A2, ,Ci, C2, as defined in (7.4), ( [7. 



( [7.5| j, ( |7. Tp anc? as acting on the function Qz{si, . . . , s^; . . . , x^^) admit 
the following expansions in z ~ 0, in terms of the operators Xk, Tk, X^i, 



4(Sm "Sm— 1) t 



-^2 = — : -i 3i + -(A-a + 2r„ + (3„_i-3s„)r_i)+0(l) 



Proof: These formulae are an immediate consequence of formulae fl7.4p , (17. 6p , 
(17.51) ■ (17. 7p . combined with the series in Lemma 7.1. ■ 
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Lemma 7.3 In the new coordinates Si, . . . , Sm, ■ ■ ■ , x^"^\ the quantities 
appearing in the basic equation liU.!^) for P„ admit the following expansions 
m z ~ 0, setting k = n/2 = 



AlCid 



J 



Im 



32a/ (s2 - Si){Sm - Sm~l) 

{ -i| - ^(A-^^Q, + 8(Si + S2- Sm-l - s^)) 

+ 2((s2-Si)2 + ( )2 + 2(si + S2)( ■^m— 1 



M7- 



m— 1 



- Si - S2)X.i + 2i'_i)T_iQ, + 0(1) 



128(S2 - Si)y/Srn " S^.i 

+ ^ ((si + S2 - - A-^i) X'_,Q, I 

+^((s™-i + s^-2(si + S2))A'_i+ 4;e_i)A'_ir_iQ, + O (^) J 

(7.9) 



512(S2 - Si)3/2^(Sm - Sm_i) 

1 



(7.10) 



256(Sm - Sm-l)A/(s2 - Si) 



+ ^A-^iT^i^ 



«4 



(7.11) 
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HI 



A2C1 + eJ 
1 



Im 



d_ 
da 



mm A 



(7.12) 



16(Sm - Sm_i)v/S2 - Si 



+ ^ (5(S2 - Sif + 24Sm_iSm + 4(S1 + S2)(5Sm_i - 3Sm)) 



2z5 



^(£i+£2) _ g^_^ _|_ 3s^^ A'_i71i 



(7.13) 



1 



64(S2 - Si){Sm - Sm-l) 



{_J_y2 (Tft _L J_ I ^X^iX^iT^i + 2X^i% + XqX^ \ 



(7.14) 



16(S2 — Si) V ~ Sm-l 

+ ^(-2X^1 + (si + S2 - (--^/--) )A'j,)Q, 
< 1 / A'-i — Xf)X-i — 2X^1% + 3X^iT^i \ „ > 

V + (1(^—1 + ^-) - 51 - S2) X_,T_J 
[ +0 (^) J 



128(S2 — Si)(Sm — Sm-l) 

' ^^^iQ. - ^Xl,T_,Q, 

+ jr{X_i + (si + S2 - Sm-i - Sm)X_^)Xl^Q, 
* ''i'js{~X^iX() — 2X^i1q + 2X-iX-iTLi > 
+ 4ri\ + 3(s^ + s^_i - si - S2)A'2,r_i)Q^ 

I +0 {^) J 
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1 



(7.15) 
(7.16) 



+ ^(-28A'_i + (18(si + S2- s^-i) + Usm)X_,)X\<i 
-16(To-l + iA'o)A'2iQ, 



+32(A'_iA'_i - r^i)r_i^ 



Proof: This is done by straightforward computation, using Lemma 7.2 and 
the series for the constants appearing in F^, H^, G^, namely 

-v2 



Jml{c) 



1 



1 + y(si + S2 - Sm-l - Sm) 



— ((S2 - Si)^ + {S„i - Sm-lf + 2(si + S2)(Sm-l + Sm)) + 0(2;" 



8 

(a-'mi) 



ml 



a 



2akJimJrr 



4:Z^{Sm - Sm-l)VS2 " Si 



-5 - y(5(si + S2) - lOSm-1 + 6Sm) 

+ ^(5(S2 - Si)2 + 24s„_iSm + 4(si + S2)(5Sm-l - Ss^)) + 0(^6) 



Proof of Theorem 0.2. 

T3 



In terms of Tj's, defined by 

= F+CiG- + F-C1G+ = CiTa - n 

= {F',F^}c, 

= F+G^ + F-Gi 

= C1T2 = F+Cfp- - F-ClF+ 

= G-C1F+ + G+CiF- , 
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the fundamental equation fl0.3l) can be written 

= -Ti ■ + T3 ■ T4 = {Ts, Talc, + T2T, 

meaning that T2, T3 and T5 only are needed, which one checks to have the 
following series in z, using Lemma 7.3, 

T = ^^-1^-1^^ I o ( 1 



1 



128(S2 - - 

+ (To - 1 + |A'o)A'2,Q, t , ^ / 1 



+ 0^ 



3 



' 1282(52 - S^fl\Sm - S™-l)2^20 + ^^19 ) ' ^ ' ^ 

Then remembering from Lemma 7.2, 

Ci = f^'^-i + 0(1) 



4VS2 - Si 

one easily computes (letting T^'^'^ be the leading coefficient of Tj in (17.171) ) 



{T2,T3}c, 

1 



{to,T30}^^+0 

1 



^31 L ^ ' 3 J Ci ' ^ V ^30 



323(S2-Sl)4(s™-S^_l)3z31 
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and then using below the trivial Wronskian relation —a^h' = {a,ab}, one 
computes 

Adding these two contributions, one finds 

= -Ti-T2 + T-i-Ti 

= {T2,Ts}c,+T2T, 
1 

lb 

+ (^To - 1 + A-^iQ, - 2(A'_iA'_i - T2i)T_iQ,| + O(^) 

(7.18) 

In [21] , Tracy and Widom show that the extended kernel for the non-intersecting 
Brownian motions tend to the extended Pearcey kernel uniformly in each 
bounded interval. Since is the log of the Fredholm determinant for that 
kernel, it follows that 

limQ, = Qo. 

Then taking the limit when 2; — >^ in (17. 18\\ leads to the PDE for Qo, ending 
the proof of Theorem 0.2. ■ 



50 



8 Appendix: evaluation of the integral over 
the full range 

Lemma 8.1 The following integral can be evaluated explicitly: 
1 



'Tkik2\ 



1=1 



i=l 



i=ki+i 

( 0<i<fei-l \ 

< 7 < n - 1 



1 .,W2^V-'"-i ciyf^yp^'^-ocvP 



= det 



v 



k\^i,A~^)) < i < fc2 - 1 / 

0<j<n-l / 



^k\ki,^ 



k\k2 



'm-\ 

n 



,-^(fcl+fc2)Jr, 



'(■^Im) 



i(fcl+fc2)2 



(8.1) 



where 



Therefore 



d 



da 



ki-l k2-l 




9a 



q; 



J] 



Im 



- (C2 + 2q; Ji^Ci) Ai log 

'^kik2 



1 a a 

:Ci- 



2 (9ci da 



kik2loga - —{ki + k2)J„ 
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Proof: From the explicit evaluation of the zero momenl0 

(27r)™/2 



V^det J- 



[Ci, . . . ,Cm-l) 

one deduces the other moments by derivation, 



/iij(±a) 



7=0 



where 



d 

A := /xoo(±a,7)"^^/ioo(±a,7) 

A_ T 

oa 



7=0 



7=0 



^4 



9 

/Uoo(±a, 7)~"^^/^oo(±a, 7) 
d 

7^ ^ii7 =F aJim 

(77 



7=0 



:= Biil) 

The following holds: 

P2j(tt) = even polynomial, P2i+i(a) = odd polynomial of a, 
which is used in equality = below, and 



^ Using 



-kiQ-,^)+{i^^)dx, ...dxm = I^^e^^O"^^'^^, 



for Q := - J-i and £ := (7, 0, . . . , 0, a). 
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where pj!*^ := and where Pi(a) are polynomials in a, independent 

of k; this feature is used in equality = below. The equality *= hinges on the 
identity 

/ i<i<fci \ 

1 < 3 < n ' 



det 



v 



[a 



1 < I < fc2 
1 < j < n 



for a constant c'^^^^ depending on ki and /c2 only Then, setting 

(27r)"^/2 



II := /xoo(±a,0) 



•v/det(-J-i(ci,...,c^_i)) 



g 2 mm CK^ 



we compute: 



U=s=u=/3=0 



= /i" det 



(^'^^•(0;)) 0<i<fei-l 
< j < n - 1 



{{-Aypj{-a)) o<i<k2-i 

< j < n - 1 



fc2(fc2-l) , ^ 

= (-1) 2 det 



/ (AVj(q;)) o<i<fci-i \ 

< J < n - 1 



{A''pj{-a)) o<i<k2-i j 



fc2(fc2-l) 

/x" (-1) 2 det 



< J < 71 

0<j<n-l 



(Pf (-«)) 





< i < fc2 - 1 / 
0<j<n-l / 



CfcifeaA*" det 



1 < 3 < n 



(((Jima?-')®) o<.<.. 



n(n— 1) 

Cfcife2/^"(Jim)^~det 



1 < 9 < 71 / 



1 < J < n / 
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(r) \nm/2 I 

A(ci, . . . ,Cm-i) 

/ m— 1 



2 



'Im 



1 



In order to evaluate the integer c^^^^' suffices to notice that this constant 
is independent of m, so that we may choose m = 1, which was done in ([5]). 
This ends the proof of Lemma 8.1. ■ 
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